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I> \ Abstract. It is showed that potentials derived from classical solutions of the SU(3) Yang-Mills 

theory can provide confining potentials that reproduce the heavy quarkonium spectrum within the 
^ | ' same level of precision as the Cornell potential. 
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In order to solve the classical Yang-Mills equations of motion, usually one writes an 
ansatz that simplifies the Euler-Lagrange equations and, hopefully, includes the relevant 
^ — dynamical degrees of freedom. In Jl|] it was proposed a generalized Cho-Faddeev- 

Niemi-Shabanov ansatz for the gluon field, where the gluon is given in terms of two 
vector fields, and Y„, and a co variant constant real scalar field n a , 
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O \ A °n= + ^fabcn%n c + Y£ (1) 

^ : 8 
Ph! with the constraints 

P,; D^n a = 0, = 0. (2) 

In it was showed that the above decomposition of the gluon field is gauge invariant 
but not necessarily complete. In the weak coupling limit, g — > 0, a finite gluon field 
requires either n = or d^n = 0. If n = 0, the gluon field is reduced to a vector field in 
the adjoint representation of SU(3) gauge group. For the other case, d^n = 0, the gluon 
is writen in terms of the vector fields A^ and Y% and includes the previous solution 
as a particular case. Accordingly, a field such that H^Oor d^n ^ does not produce 
a finite gluon field in the weak coupling limit and, in this sense, can be viewed as a 
nonperturbative field. Among this class of fields, the simplest parametrisation for the 
covariant scalar field 1 is rf = 8 ai {— sin0) + 5 a2 (cos0). Then 

A c l = n a A^ + 8 a3 ^0 + 5 fl8 C M , (3) 

where = Y^. The classical Lagrangian and equations of motion are independent of 
and are abelian like in A^ and C^. Among the possible nonperturbative gluons given 



1 From the constraint equation Dn = it follows that n 2 is constant. Our choice was n 2 = 1 . A different 
value for the norm of n is equivalent to a rescaling of A. 



by ([T]), the simplest configuration has A = C = 0. The coupling to the fermionic fields 
requires only the Gell-Mann matrix A 3 , decoupling the different colour components. 
This suggests, naively, that such a field is able to produce either confining, non-confining 
or free particle solutions for the quarks. 

The classical equations of motion are independent of 6. However, a choice of a gauge 
condition, provides an equation for this field. For the Landau gauge, 6 verifies a Klein- 
Gordon equation for a massless scalar field. Note that there is no boundary condition 
for 0, i.e. the usual free particle solutions of the Klein-Gordon equation are not the only 
possible ones. Indeed, writing 9(t,r) = T(t)V(f), then 

T"{t) _ VV(F) _ A 2 (A , 
Tit) ~ - A > U , (4 ) 

T(t) = ae M + be~ M , (5) 
V (f) = E Vi{r) Y lm (Q) , (6) 

Vi{r) = f z h + i,i{z) + ^K l+l/2 (z) , (7) 

where z = Ar and //+i/2(z) and ^+1/2(2) are modified spherical Bessel functions of the 
I s *- and l rc ^ kind 2 . The lowest multipole solution is 

. , sinh(Ar) e~ Ar 

V (r) = A y —L + B (8) 

r r 

and the associated gluon field is given by 

Aq = A{e At -be At ) V {r), (9) 
Aq = -{e^-be**) Wb(r). (10) 

From the lowest multipole solution one can derive a potential, which maybe suitable 
to describe heavy quarkonium. Indeed, assuming that quarks do not exchange energy, in 
the nonrelativistic approximation and leading order in 1/m, the spatial function in Aq, 
Vb(r), can be viewed as a nonrelativistic potential 3 and one can try to solve the associated 
Schrodinger equation. For the potential (HJ), the wave function goes to zero faster than 
an exponential for large quark distances, 

yr(r) = exp I J— exp ( j }. (11) 
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As a first try to compute the heavy quarkonium spectra, we fixed A, B and A minimis- 
ing the square of the difference between Vo(r) + Constant and the Cornell potential |0] 
Vcornell = e/r + or (e = —0.25, \J~G = 421 MeV) integrated between 0.2 fm and 1 fm. 



Note that, by definition, the mass scale A is independent of a rescaling of the gluon field. 
The potential is ~ l/r for short distances and goes to infinity for large quark distances. 
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FIGURE 1. Heavy quarkonia spectra in MeV. The plots include the spin averaged experimental values. 



This optimisation provides the following parameters A = 5.4, B = — 1 .0, A = 28 1 MeV, 
Constant = - 1 190 MeV; for these values -2AM eV < V Cor neii ~ (Vb + Constant) < 64 
MeV in the integration range considered. Then, we can compare the Schrodinger equa- 
tion spectrum for the charmonium (m c = 1 .25 GeV) and for the bottomonium (m^, = 4.25 
GeV) for the two potentials. The spectrum for the new potential shows an equal level 
spacing for both the charmonium and bottomonium spectra. If the Vq charmonium spec- 
trum is quite close to the Cornell spectrum, the botomonium shows clear deviations; 
see figure [T] The differences are the result of overestimating the strengh of Vo( r ) for 
smaller distances. Indeed, one can improve our potential linearising the full QCD equa- 
tions around the above configuration. To lowest order, this is equivalent to add a term 
like k/r to Vq. Computing k perturbatively 4 adjusting the M[(1P)] —M[(\S)] bottomo- 
nium mass difference, gives k = 0.2448251. The heavy quarkonia spectra, including this 
correction, is given in figure [T] 

In conclusion, classical configurations seem to be able to produce a spectra close to 
the Cornell potential. Hopefully, this is an indication that these configurations can be of 
help to understand strong interaction physics. Of course, there are a number of issues that 
need to be further investigated (definition of the potential parameters, inclusion of time 
dependence, decay rates). We are currently working on these topics and will provide a 
report soon. 
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For each stationary, the shift in energy due to this term is compatible with a perturbative treatment. 
Corrections are clearly below 10-20%. 



